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Abstract 

The cluster reduction method for the Yakubovsky equations in configu¬ 
ration space is used for calculations of zero-energy scattering in four-nucleon 
system. The main idea of the method consists in making use of expansions for 
the Yakubovsky amplitudes onto the basis of the Faddeev components for the 
two-cluster sub-Hamiltonian eigenfunctions. The expansions reduce the origi¬ 
nal equations to ones for the functions depending on the relative coordinates 
between the clusters. On the basis of the resulting equations the N—(NNN) 
zero-energy scattering problems are solved numerically with the MT I III model 
for N-N forces and neglecting the Coulomb interaction between protons. 


1 Introduction 

In this paper we continue our investigations of the four particle scattering prob¬ 
lem making use the cluster reduction method for Yakubovsky equations in con¬ 
figuration space formulated in [1], ||]. The Yakubovsky differential equations 
(YDE) being the direct generalization of the three body Faddeev approach, 
inherit all advantages of latter. Two main of them are the very simple way to 
impose the boundary conditions corresponding to all elastic and rearrangement 
channels and minimum inputs needed to formulate the equations and boundary 
conditions. 

For the first time, YDE have been used for calculations of the four identical 
particles bound state in Q . It was observed that even for the simple model N — 
N interactions the converged results within the usual numerical approximations 
of the YDE could be obtained with supercomputer facilities only. So, 
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the progress in four nucleon calculations with YDE depends on a new numer¬ 
ical solution technique. For the bound state problem such a technique was 
elaborated in || and allowed to improve the results of ||] considerably. For the 
scattering problem one of the new solution method was introduced in |l], |], . 

Here, we employ the method for calculations of zero-energy N-(NNN) 
scattering neglecting the Coulomb interaction between protons. This system 
could be considered as a model for n- 3 He, p- 3 H or n- 3 H scattering depending 
on a value of the total isospin. 


2 Yakubovsky equations 

In the Yakubovsky approach, the four particle wave function T should be de¬ 
composed into components in one to one correspondence to all chains of par¬ 
titions. The chains consist of two cluster 02 ( e.g ., ( ijk)l or ( ij)(kl )) and three 
cluster 03 (e.g., (ij)kl) partitions obeying the relation 03 € ci 2 - The latter 
means that the partition 03 can be obtained from partition ci 2 by splitting of 
one subsystem. It is easy to see that there exist 18 chains of partitions for the 
four particle system. 

The Yakubovsky wave function components can be defined by formulas [|j, £|] 
^0,3(12 ^ 0,3 (E + )V a3 Y Ro(E + )V h V, ( 1 ) 

(b 3 ^Q3)Sa 2 


where 


Ra 3 (E) = (Ho + V a3 - E)-\ Ro(E) = (H 0 - E)~\ E + = E + iO, 

and V as stands for the two particle potential acting inside the two particle 
subsystem of a partition 02 - If the function is the solution of the Schroedinger 
equation 

(H 0 + Y / Va 3 -E)* = 0, 

a3 

then the components \k a3 a 2 obey the Yakubovsky equations |3], 0 
(Ho T Va 3 E )* asa2 + V a3 'y ( * C3a2 = 

(c3^Q3)Sa2 

-K 3 E £ ***■ ( 2 ) 

d,2^a,2 (d 3 ^a 3 )£a2 

There exists a remarkable rule which allows to construct the Schroedinger 
equation solution from the components ^a 3 a 2 §0 

* = ££*««■ 

«2 a 3 

Note that, the only inputs needed to formulate the YDE (|2|) are the inter¬ 
particle potentials. 
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The second advantage of the YDE consists in the structure of the bound¬ 
ary conditions. In contrast to the total wave function for the multichannel 
scattering, the asymptotics of the Yakubovsky component has the one channel 
form. Namely, T asa2 involves in the asymptotical region the characteristics of 
the bound states for subsystems of the partitions <22 only. As a consequence, 
the asymptotic behavior of T asa2 can be described in terms of one set of the 
relative coordinates corresponding to the chain 0302- 

Let us turn to description of this asymptotics. To this end introduce the 
relative Jacobi coordinates for the four particle system. There exist 18 sets 
of Jacobi coordinates in one to one correspondence to 18 chains of partitions. 
Note, that among them only two sets are topologically different. One corre¬ 
sponds to chains with partitions 02 of the type 3 + 1 (e.g., (123)4) and the other 
one does to chains with partitions <22 of the type 2 + 2 (e.g., (12)(34)). These 
coordinates can be given by formulas 

x „ 3 = r i- rj 

y a 3 a 2 = (r i + rj)/2 - r k (3) 

Za 2 = (r* + I -j + r fc )/3 - 17 

for 02 = (ijk)l, 03 = (ij)kl and 

x fe3 = r, - r j 

yb 3 b 2 = r k - r 1 (4) 

z b 2 = ( r i + O')/ 2 - ( r fc + r l)/ 2 

for b- 2 = ( ij)(kl),b 3 = (ij)kl. 

In this paper we are considering only binary processes, it assumes that only 
the two clusters channels are open. In this case, the components which corre¬ 
spond to the initial state with binding clusters of subsystems of the partition 
I 2 have the form 


^a 3 o 2 (X, P/ 2 ) = da 2 l 2 ^al(^a 2 )eX.p{i(p a2 ,Z a2 )} 

+Val( X tt2 )Ua 2 l 2 ( z a2 5 P/ 2 ) J - ^^a 2 a3,/ 2 (X, p/ 2 ). 

Here x a2 = {x a 3 ,y a 3 ,a 2 } 5 X = {x a 2 ,z a2 } and p ^ 2 is the relative momentum 
conjugated to the vector z i 2 . The momentum p^ 2 and the energy of the system 
E obey the equation 

n 2 

E £+ + (P^ 2 ) • 

m 

s a2 and are binding energies and the Faddeev components of the bound 
states for the two cluster Hamiltonians 

hci 2 = h® 2 + ^2 Va 3 = — — (^x a3 + ^a 2 Ay a3 a 2 ) + ^ K 3 ( x a 3 ) 5 

a 3 £a 2 a 3 £a 2 

where d 02 = 3/4 for the partitions 02 of the type (3+1) and d a2 = 1 for the 
partitions 02 of (2+2) type. 


3 



(6) 


Amplitudes U become the spherical waves 

TT , x „ ex.p{iy/E - e a 2 \z a2 \} 

U (I2I2 (^0,2 i Pi2) ~ -Aa2(2 I I (. 1^021 * OOj. 

z a 2 I 

Amplitudes SU a2 .a 3 ,i 2 are exponentially decreasing functions as |X| —> 00 cor¬ 
responding to virtual breakup processes. 


3 Cluster reduction of YDE. 

The cluster reduction of YDE j|] consists in expansion of the components T a3 a 2 
onto the basis consisting of functions depending on the coordinates x a2 . The 
most suitable basis for that is the one of the eigenfunctions for the operators 
from the left hand side of the YDE (|2|) 

K + v.,)€h + v m £ P) 

(C3^a 3 )£a2 

The expansion of the components has the form 

00 

*a 3 a 2 = T,€lk(Xa 2 )F* 2 M, ( 8 ) 

k=0 

where the coefficients ( z a2 ) depend on the relative coordinates between the 
clusters only. 

The basis of the solutions of the eqs. (|7|) is not an orthogonal one, be¬ 
cause the eqs. © are not Hermitean. Fortunately, it is easy to show that the 
Hermitean conjugated equations to the eqs. are of the form || 

« + *WO + £ (f) 

(C 3 ^“ 3)£“2 

In the papers || it is shown that equations ([?]) and (0) allow the solutions of 
two classes. They are the physical ones, i.e., the solutions of eqs. © which 
sum 

Y VC,* = ^2 ,k 

036 a 2 

is the solution of the Schroedinger equation 

ha,2'4 ) a2,k ^ a2^ a 2M’’ 

and the solutions of eqs. (p) with the property 

C* = 

for any a 3 G ai- The second class of solutions of eqs. is formed by the 
spurious solutions for which 

£ O = o- 

I3& a 2 
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Unfortunately, there is no a such simple formula for the spurious solutions of 
eqs. ([)[). It is important to note, that the complete basis is formed only by the 
physical and spurious solutions together. Complete bases of solutions of eqs. 
® and (D are biorthogonal 

E = «»• 

I3& a 2 

This property allows to apply the projection technique in order to obtain the 
equations for coefficients F* 2 from series (||) after substitution of eqs. (j|) into 
eqs. (||). The result reads 

(~S^Az, 2 -E + e* 2 )F« 2 (z a2 ) = 

- E E E E 4%AM), (10) 

a 3Ga2 d 2 ^a 2 (rf3^a3)ea2 />0 

where the brackets (.|.) stand for the integration over x a2 . 

The asymptotic boundary conditions for the coefficients F a 2 M can be 
obtained from (jsj) , (j^) by projecting and have the following ’two body’ form 

F a 2 ( z a 2 ) ~ <W 2 expi(p a2 ,z a2 ) + ^ a2 ; 2 |z a2 | _1 expiy / E - e a2 |z a2 | 
for the open channels and 

F a 2 ( z a 2 ) ~ 0 ,k> 1 


for the closed channels. 


The equations (10) are the desired coupled channel equations for the two 
cluster collisions in the four body system. 


4 Application to the N — NNN scattering 
problem. 

The eqs. 0 ) can be applied to solving the IV — NNN scattering problem after 
a suitable partial wave analysis. We have used the MTI-III model for the N—N 
forces and the isotopic spin formalism. Within this conditions and neglecting 
the Coulomb interaction the channels corresponding to total value of isospin T 
are independent. 

The calculated values of scattering lengths in N — NNN system for a given 
values of the total spin S and the total isospin T are given in the Tab. 1 with 
data of direct solution of YDE from the paper [(J. As one can see the agreement 
of calculations is rather reasonable. In the Tab. 2 we collect the available data 
for n — 3 H (T = 1) scattering lengths calculated within various approaches. As 
one can see the agreement of the calculations is quite good. By the end in the 
Tab. 3, we give the results of calculations of scattering lengths in the channel 
T = 0 performed under conditions when the spurious part of basis functions is 
removed from expansions (pj) in order to show how the presence of the spurious 
solutions affects on the results of calculations. 
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Tab. 1 


s 

T 

[ 3 ] 

present work 

0 

0 

14.75 

14.7 

1 

0 

3.25 

2.9 

0 

1 

4.13 

4.0 

1 

1 

3.73 

3.6 
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Tab. 2 


Refs. 

^S = o fm 

As=i fnr 

Present work 

4.0 

3.6 


Too 

4.09 

3.61 


nj 

4.23 

3.46 


0 

3.905 

3.597 


Hi 

3.8 

4.9 

1 

n 

3.38 

3.25 

m 

4.13 

3.73 

exp. 15] 

3.91 ± 0.12 

3.6 ± 0.1 


















Tab. 3 


□ 

A 

3+1 phys. 

3+1 spur. 

2 + 2 phys. 

2 + 2 spur 


14.7 

+ 

+ 

+ 

+ 

Q 

15.4 

+ 

- 

+ 

- 

i 

2.9 

+ 

+ 

+ 

+ 

i 

1.4 

+ 

+ 

+ 

- 

i 

1.2 

+ 

- 

+ 

- 
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